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In this paper, we investigate topological properties of the ground state of the SU(A'^) Heisenberg chain, which 
is argued to be relevant to the Mott-insulating phase of alkaline-earth cold fermions in a one-dimensional optical 
lattice. By calculating the entanglement spectrum, we show that the ground state is in one of the topological 
phases protected by SU(A'^) symmetry. We then discuss an alternative characterization of it with non-local string 
order parameters. We also consider how the reduction of the protecting symmetry affects the topological phase 
paying particular attention to the entanglement spectrum. 
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I. INTRODUCTION 

Symmety in physics not only is the key to understand¬ 
ing phases of matter but also play a vital role in unifying 
seemingly different things and uncovering fundamental prin¬ 
ciples underlying them. In particular, unitary groups have 
been playing very important roles in quantum mechanics as 
the orthogonal groups in classical mechanics. For instance, 
SU(3) is the fundamental symmetry underlying the quan¬ 
tum chromodynamics (QCD) of strong interactions. In tra¬ 
ditional condensed-matter physics, however, high symmetry 
like SU(iV) is usually realized, aside from few exceptions, 
only in rather idealized situations and has been mainly used as 
mathematical convenience that makes problems tractable. For 
instance, in the large-iV approximations, we replace the phys¬ 
ical symmetry SU(2) with SU(7V) and use 1/7V as the (small) 
control parameter of the approximation hoping that there is a 
smooth crossover down to N = 2. 

Recent suggestions*'^ that SU(iV)-symmetric fermion sys¬ 
tems could be simulated using the alkaline-earth atoms and 
their cousins (*''*Yb, ®^Sr, etc.) loaded in optical lat¬ 

tices opened a new era of SU(iV) physics^’"* (see, e.g.. Refs. 5 
and 6 for recent reviews). For instance, the S\J(N) gener¬ 
alization of quantum magnetism is of direct relevance to the 
Mott-insulating regime of these systems. The SU(7V) “spin” 
models provide us with examples of underconstrained systems 
that yield, on top of usual “magnetically ordered” states, var¬ 
ious unconventional states, e.g., deconfined criticalities,^’^ an 
algebraic spin liquid® and a chiral spin liquid.*** 

On the other hand, topological states of matter* * have been 
subjects of extensive research for the past decade. Since the 
advent of topological insulators and superconductors,*^ it has 
been widely realized that there exists a special class of “topo¬ 
logical” phases that is stable only in the presence of certain 
symmetries.*** *^ This class of topological phases is called 
“symmetry-protected topological (SPT)”*^ as it is topologi¬ 
cally protected only when we impose symmetries on the sys¬ 
tem in question, and otherwise they reduce to trivial ones. The 
catalogue of possible topological phases depends crucially on 
the symmetry we impose and different lists of possible phases 


may be obtained for different protecting symmetries (see, e.g.. 
Ref 18 for a catalogue of SPT phases). One defining property 
of SPT phases is the existence of gapless boundary excitations 
(edge states) that are intrinsically different from those in the 
gapped bulk. A modern mathematical way of observing the 
edge states would be to use the entanglement spectrum*® that 
is obtained solely from the ground-state wave function. In the 
following, we heavily use the entanglement spectrum in char¬ 
acterizing topological phases. 

Despite the recent effort*^ in systematically enumerating 
possible SPT phases in one dimension, not much is known, 
except for a few examples, about how to observe those phases 
in realistic settings. Recently, it has been suggested^**’^* that 
a class of SPT phases is realized in the Mott-insulating re¬ 
gion of the alkaline-earth cold fermions, and this is one of 
the motivations of our study here. Specifically, deep inside 
the Mott phase at half-filling, the low-energy physics of a 
system of alkaline-earth fermions is described by an S\J{N) 
“spin” model (see Secs. IIA and IIB) whose ground state 
is expected to be in one of the topological phases predicted 
in Ref. 22. Therefore the alkaline-earth fermions provide us 
with a unique arena for the realization of new SPT phases in 
a very controlled manner. Our goal is to clarify the nature 
of the ground state of the above S\J{N) spin Hamiltonian in 
several complementary ways and demonstrate the use of non¬ 
local string order parameters to detect the phase. 

The outline of this paper is as follows. In Sec. II, we 
introduce the S\J(N) Heisenberg model and sketch how 
it is derived as the effective Hamiltonian for the Mott- 
insulating phase of the alkaline-earth cold fermions on a one¬ 
dimensional optical lattice. A variant of the Heisenberg model 
that gives useful insights about the topological properties of 
the original model is introduced as well. After briefly sum¬ 
marizing the minimal background of SPT phases expected for 
our SU(A^) spin systems, we try, in Sec. Ill, to characterize 
the topological properties of the ground state of the SU(A^) 
Heisenberg model using its entanglement spectrum. By care¬ 
fully investigating the structure of the spectrum obtained for 
N = 4, we present a strong evidence that the ground state 
of the SU(4) Heisenberg model is in one of the SU(4) topo- 
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logical phases. In Sec. IV, we present an alternative way of 
characterizing the SU(7V) SPT phases using non-local string 
order parameters. 

Although the alkaline-earth fermions, that motivated our 
study, possess very precise S\J{N) symmetry, it would be 
interesting theoretically to consider the situations where the 
original S\J{N) symmetry gets lowered. We investigate this 
problem in Sec. V to find that, depending on N, the system re¬ 
mains topological even after the SU(A^) symmetry is relaxed. 
Summary of the main results is given in Sec. VI. 


II. MODEL 


In this paper, we consider the ground-state properties of the 
following Hamiltonian 

N^-l 

^Heis = J E 
A^l 


where Sf (A = 1,..., — 1) denote the S\J{N) genera¬ 

tors. In S\J(N), instead of fixing spin S, one has to specify 
the irreducible representation(s) to which the generators 
belong. In the following, Sf ( A = 1,..., — 1) denote, un¬ 

less otherwise stated, the SU(A^) generators in the irreducible 
representation characterized by the following Young diagram 
with N/2 rows and two columns: 


N/2 



{N = even) . 


( 2 ) 


It is well-known that the low-energy physics of the SU(A^) 
Heisenberg model depends crucially on the representation(s) 
we put on the individual lattice sites. For the fully- 
symmetrized representation • • • □ (ric boxes), the exact 
Bethe-ansatz solutions are available;^^“^^ the ground state is 
known to be gapless and described by the level-ric SU(A^) 
Wess-Zumino-Witten conformal field theory with the central 
chargee = nc(W^ — l)/(W-|-nc)-^® For any translationally in¬ 
variant choice of representations (i.e., the same representation 
is assigned on every site), we can show that the S\J{N) chain, 
which has a unique (finite-size) ground state^^, is either gap¬ 
less or has degenerate ground states (with broken symmetries) 
provided that the number of boxes ny in the Young diagram 
is not divisible by In other words, except for the cases of 
ny = 0 (mod N) [including the one shown in Eq. (2) which is 
relevant to our spin chain], this statement excludes the possi¬ 
bility of gapped topological ground states. Remarkably, this is 
perfectly consistent with the recent group-cohomology clas¬ 
sification of the gapped SPT phases^^ (see Sec. IIIB for the 
detail). There is also an attempt^^ at summarizing these ob¬ 
servations into a “generalized” Haldane conjecture. 

Some insights about the nature of the ground state of (1) 
are gained from the large-analysis^°“^^ as well. For N/2 
rows but with a single column, the ground state is expected 
dimerized,^® while, for two columns, we may have a gapped 
translationally invariant ground state,which we will argue 
to be topological. 


A. Relation to cold fermion systems 


It has been argued in Refs. 20 and 21 that the Hamilto¬ 
nian (1) emerges as the effective Hamiltonian in the Mott- 
insulating region of the alkaline-earth cold fermions loaded 
in a one-dimensional optical lattice at half-filling. To empha¬ 
size the relevance of our results to experimentally realizable 
systems, we sketch how the model ’Hneis is derived from the 
cold-fermion systems in the Mott region. 

It is known that the decoupling between the nuclear spin 
(/) and the total electron angular momentum makes it possi¬ 
ble to organize the {21 -f 1) nuclear-spin states of each atom 
into a multiplet of larger SU(2/ -f l)-symmetry. Specifically, 
the interaction between two like alkaline-earth atoms does 
not depend on the nuclear-spin states of each and hence is 
SU(2/ -f l)-symmetric.'’^ Moreover, one can add one more 
degree of freedom (orbital) by taking into account the first 
meta-stable excited states (in ^Pq', denoted as “e”) as well as 
the atomic ground state in ^S'o (“g”)/^ That this SU(2/ -f 1)- 
symmetry holds for both orbitals with very high accuracy has 
been verified in recent scattering-length measurements. 

When loaded into a one-dimensional optical lattice, the sys¬ 
tem of alkaline-earth cold fermions is described by the follow¬ 
ing Hubbard-like Hamiltonian^ 


N 

'^G=-Y1 E E +h.C.) 


i m—g,e a—1 


u, 


(m) 


- E Fg E’^^’^+E E -A^'nni,t{nra,i - ^) 

m—g,e i 




I m—g,e 

‘^ga, N'eP, Ng0, Nea,ii 


l,OLp 


(3) 


where N = 2/ -f 1 denotes the number of nuclear-spin 
states and the operator ^ creates an atom in the internal 
state (a, to) (a = 1,..., N, m = g,e) at the site i. The 
number operators are defined as nma,i = Nma,i and 
i = 'l2a=i As the two orbitals are not symmetry- 

related, the hopping amplitudes (to = g, e), the chem¬ 
ical potential and the intra-orbital interaction in 

general are different for the two orbitals. The inter-orbital ex¬ 
change (or, Hund coupling) is crucial in determining the 
nature of the Mott-insulating phases.^' 

Clearly, the Hamiltonian (3) is invariant under the SU(A^) 
transformation 


N 

Gma,i t E [Z^GSU(7V)] (4) 

/3=1 

as well as the multiplication of a global U(l) phase: 

Gma,i t Q C^q, ^ . (5) 

Borrowing a terminology from the electron systems, we call, 
in the rest of this paper, the degree of freedom associated 
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with (5) “charge”, although the fermions Cma,i are charge- 
neutral in the cold-atom context. This and the related systems 
have been investigated extensively both for SU( 2 )^®“^® and for 
SU(iV).20-21.40 


irreducible representation specified by the Young diagram (2). 
Therefore, one sees that the model "HHeis [eq.(l)] describes 
the low-energy physics of the alkaline-earth cold fermions 
[eq.(3)] in the Mott-insulating phase (for J = 14® ® > 0). 


B. Strong-coupling limit 


C. Solvable Hamiltonian 


Recently, it has been argued^®’^' that for large positive (7^"*^ 
and there exists a topological Mott phase protected by 
SU(4)-symmetry."^* In order to consider the Mott-insulating 
phases, it is convenient to start from the strong-coupling limit 
U^\ Vg, ^ In this limit, charge fluctuations are 
strongly suppressed and the SU(iV) “spin” and orbital domi¬ 
nate the low-energy physics. One may introduce the psuedo- 
spin operator 7;“ = 5 Ea,/ 3 .m (a = x,y,z) 

for each orbital to rewrite the single-site (i.e., = 0) part 

of the Hamiltonian as 

= 0 ) = ha,om.c(z) 
i 

^atomic(^) = ^ (Mc “t” f-^g) “f ( 6 ) 

— {y,g — fie) + UdiffT^ni 

with the following coupling constants 

u = + 2^0), C/diff = - ut), 

J = Vr, J. = 1{U^^^ + U^^^-2 Vg), (7) 

Mm = ^ (2 Mg"^ + ^4™^ + Kx {'m = g,e) . 

Let us consider the case of half-filling where each site is 
occupied by N fermions on average. The Fermi statistics al¬ 
lows {2N)\/{Nl)'^ states and, out of them, the optimal ones 
are chosen by the orbital-dependent terms [the last four terms 
in /Ltomic(*)]; when N is even and I4f’® is positive, the states 
that transform under SU(A^) as the irreducible representation 
(2) are the ground states of When iV = 4, they 

form the 20-dimensional representation of SU(4). For these 
states, the orbital pseudo-spin is quenched and only the 
SU(iV) degree of freedom remains. When N is odd, on the 
other hand, both SU(iV) spin and orbital are active and we 
obtain, in general, SU(A^)-orbital-coupled models. In the fol¬ 
lowing, we consider only the case with even-7V where pure 
spin models are obtained. 

Interactions among the remaining SU(iV) spins are derived 
by the second-order perturbation in as^* 

2 \ t/ -f (7d.ff + J + i [/ - C/d.ff + J + i J * ’ 

( 8 ) 

where we have introduced a short-hand notation Si ■ Si+i = 
being the SU(A^) generators in the 


Unfortunately, the Heisenberg Hamiltonian (1) cannot be 
solved exactly. However, one can design a solvable model 
Hamiltonian whose ground state may share important proper¬ 
ties with that of the original Heisenberg model (1). Clearly, 
when N — 2, the Affleck-Kennedy-Lieb-Tasaki (AKLT) 
model proposed in Refs. 42 and 43 will do the job; 


n 


N=2 

VBS 



1 

3 


(S. • S,+i)' 


(9) 


where denote the spin-1 operators. Its (rigorous) 
ground state, dubbed the valence-bond solid (VBS) state, is 
constructed'^^’"^^ by hrst decomposing an S' = 1 on each site 
into a pair of S = l/2s, forming uniform tiling of dimer sin¬ 
glets (‘valence-bond solid’) among the neighboring sites, and 
then fusing the S = 1/2 pairs back to the original spin-ls. 

Suggested by the above construction of the VBS ground 
state, we can think of constructing the model ground state by 
first preparing two auxiliary ‘spins’ 


N/2 



{N = even) 


( 10 ) 


on each site and pairing such spins on the adjacent sites into 
SU(W) singlets (see Fig. 1). The VBS ground state is ob¬ 
tained by projecting the product of the two fictitious spins 
on each site onto the physical Hilbert space characterized by 
the Young diagram in (2) (see Fig. 1). In the following, we 
call this kind of states the SU(A^) VBS states.The parent 
Hamiltonians for these states read, e.g., for N = ^nd 
for W = 6 as"^^ 


nj N—4. 

^VBS 

= ^ • 5i+i -I- {Si ■ -I -{Si ■ Si+iY 


( 11 ) 


and 


‘^*+1 + 


508 




( 12 ) 


respectively.^^® (In writing down the above expressions, we 
have normalized the generators Si in such a way that the 
lengths of the simple roots are all v/2-) The dimensions of 
the physical SU(W) ‘spin’ multiplet on each site are 20 and 
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175 for = 4 and 6, respectively. In Refs. 20 and 21, the 
ground state wave function of has been obtained in a 

matrix-product-state (MPS) form (see Appendix A). Clearly, 
the higher-order terms are rapidly suppressed as we go to 
larger-A^. This suggests that the larger N is, the better the 
VBS state shown in Fig. 1 approximates the ground state of 
the original Heisenberg model (1). This is quite natural in 
view of the large-A^ results.These models will serve as an 
ideal starting point for the study of the topological properties. 




ffl 




(a) 


(b) 



SU(4)-singlet 




SU(N)-singlet 


>N/2 


FIG. 1. (Color online) (a) Ground state of SU(4) VBS model 
[eq.(ll)]. Two 6-dimensional representations (‘fictitious spins’) are 
projected onto a physical 20-dimensional representations, (b) Simi¬ 
lar construction applies to the cases with larger N as well. 


III. SYMMETRY-PROTECTED TOPOLOGICAL PHASES 


In this section, we try to characterize the nature of the 
ground state of the S\J{N) spin chain (1). Specifically, in 
Sec. me, we show that the ground state of the model (1) 
shares essentially the same properties with that of the solvable 
VBS models and that it is in fact in one of the SPT phases. 
Being topological, this class of topological phases defies the 
traditional characterization with broken symmetries and the 
associated local order parameters. One way is to use the physi¬ 
cal edge states to distinguish between topological phases from 
trivial ones. However, this approach is not quite satisfactory in 
the following respects. First, even topologically trivial states 
may have certain structures around the edges of the system, 
as, e.g., the spin-2 Heisenberg chain does.^^’^^ Second, in or¬ 
der to see the edge excitations, it is necessary to consider the 
excitation spectrum, while the topological properties are in¬ 
trinsic to the ground state itself and should be seen only by 
examining the ground-state wave function. 

Recently, the use of the entanglement spectrum in char¬ 
acterizing topological phases has been suggested in Ref. 19. 
This is based on the observation that the entanglement spec¬ 
trum resembles the spectrum of the physical edge excita¬ 
tions. The idea has been successfully applied to various 
systemsand enabled us to characterize topological 
phases and quantum phase transitions among them. In this 
section, we present a clear evidence from the entanglement 
spectrum that the ground state of the SU(4) Heisenberg model 


(1) is indeed in the SPT phases protected by SU(4) [PSU(4), 
precisely] symmetry. 


A. Haldane phase -an SPT primer 

To understand the nature of the SPT phases in the case 
of SU(A^) symmetry, it is convenient to begin with the sim¬ 
plest case N = 2. In 1983, Haldane conjectured^"^’^^ that 
the ground-state properties of the spin-S' Heisenberg chain 
are qualitatively different according to the parity of 25'; when 
25 = even, the ground state is in a featureless non-magnetic 
phase (Haldane phase) with the gapped triplon excitations 
in the bulk, while, for odd 25, we have a gapless (i.e., al¬ 
gebraic) ground state with spinon excitations. This conjec¬ 
ture has been later confirmed both by the construction of a 
rigorous example^^’"'^^’^® [Eq. (9)] and by extensive numerical 
simulations.Soon after, it has been pointed out that the 
featureless gapped ground state of the integer-5 spin chains 
may have a hidden “topological” order characterized by non¬ 
local order parameters®*'^^ at least when 5 is an odd integer.®^ 

However, it was not until the concept of SPT phases was 
established that the true meaning of “topological order” in 
the Haldane phase was understood.'^ Now it is realized that 
the gapped phases in integer-spin chains with some protecting 
symmetry (e.g., time-reversal, reflection) are further catego¬ 
rized into topological phases and the other trivial ones. To un¬ 
derstand the difference, it is useful to consider how the ground 
state in question transforms under the symmetry operation. As 
the ground state is assumed symmetric, the bulk does not re¬ 
spond to the symmetry operation but the edges do. As the 
consequence, the symmetry operation gets fractionalized into 
two pieces; one acts on the left edge and the other on the right. 
For instance, the VBS ground state |5 = 1 VBS)a _/3 of the 
spin-1 AKLT model (9) hosts two emergent 5 = | spins (i.e., 
a, P =t, i) on both edges and hence transforms under the 
SO(3) rotation as 

|5 = 1 VBS)„,^ ^ = 1 VBS)„.,;3. , 

a' ,(3' 

(13) 

where U is the 5 = ^ rotation matrix of SU(2). Putting it an¬ 
other way, U serves as the mathematical labeling of the phys¬ 
ical edge states. It is important to note that U in the above in 
general is a projective representation of SO(3) as both f/l and 
U appear simultaneously in the equation. 

Since this U belongs to a non-trivial projective repre¬ 
sentation that is intrinsically different from any irreducible 
representations of the original SO(3), one sees that |5 = 
1 VBS)q is in a non-trivial topological phase with emergent 
edge states. On the other hand, one can construct another ex¬ 
act ground state of a spin-1 chain which transforms as above 
but with U belonging to the spin-1 representation. Since the 
spin-1 representation is trivial in the sense of projective rep¬ 
resentation of SO(3), one can kill the would-be edge states 
by continuously deforming the Hamiltonian^® and this ground 
state is in a trivial phase. This reasoning may be readily gen¬ 
eralized; when U transforms like a half-odd-integer spin, the 
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phase is topological, while when U transforms in an integer- 
spin representation [i.e., linear representation of SO(3)], the 
system is in a trivial phase. What is crucial in the topological 
properties is not the bulk spins at the individual sites but the 
edge spins. 

For later convenience, we summarize the situation in terms 
of Young diagrams. The spin-S” representation of SU(2) is 
represented by the following Young diagram: 

m-'-D- (14) 

2S boxes 

With this in mind, the above result may be summarized as 
follows; when U belongs to the representations 

n, . (15) 

the state represented by the corresponding MPS is topologi¬ 
cally non-trivial, while the phase is trivial for U transforming 
in 


(16) 


That is, the number of boxes (mod 2) in the Young diagram for 
the representation to which U belongs labels the topological 
classes protected by SO(3) and leads to the Z2 classification 
of the SO(3)SPT phases. 

B. SU(Ar) topological phases 

Using the MPS representation®^ of the gapped ground state 
in one dimension, the above “physical” idea can be general¬ 
ized and made mathematically precise. In fact, when a given 
ground state that is represented by an MPS 

A{mi)A{m2) ■ ■ ■ A{mL)\mi)<^ ■ ■ ■ (g)|TOL) (17) 

{rrii} 

is invariant under some symmetry G, a Z7-dimensional unitary 
matrix Ug (g G G) exists such that®® 

A{m,)^e^^^UlA{m,)Ug, (18) 

where A{mi) denotes the DxD MPS matrices corresponding 
to the local physical state \mi) and is a phase that de¬ 
pends on G. As has been mentioned above, the unitary matrix 
Ug is in fact a projective representation of the symmetry G, 
that corresponds to the physical edge states.*^ Therefore, the 
enumeration of topologically stable phases in the presence of 
symmetry G boils down to counting the possible (non-trivial) 
projective representations of G.'® 

This problem was solved for SU(A^) and other Lie groups in 
Ref. 22 and the picture in the previous section basically gen¬ 
eralizes to the case of S\J{N) with some mathematical com¬ 
plications. Now the role of SO(3) in the previous section is 
played by PSU(A^) ~ SU(N)/Zn [note SO(3) ~ PSU(2)]. 
Considering PSU(A^) instead of S\J{N) amounts to restrict¬ 
ing ourselves only to the irreducible representations of S\J(N) 


specified by Young diagrams with the number of boxes ny di¬ 
visible by N [i.e., ny = Nk (fc = 0,1,...)]. This subset of 
irreducible representations roughly corresponds to the integer- 
spin ones in the SU(2) case. As in the previous section, the 
topological class of a given ground state (typically written 
as an MPS) is determined by looking at to which projective 
representation the unitary Ug of the state belongs. Since in¬ 
equivalent projective representations of PSU(A^) are labeled 
by uy (mod N),^^ there are A^ — 1 non-trivial topological 
classes (as well as one trivial one) specified by the la¬ 
bel ritop = ny (mod N). In the following, we use the name 
“class-ritop” for these topological classes (the class-0 corre¬ 
sponds to trivial phases). For instance, one can readily see 
that the “VBS states” (which are different from ours) investi¬ 
gated in Refs. 43, 67, and 68 fall into the class-1 and A^ — 1 of 
the PSU(A^) SPT phases (see Supplementary Material). Quite 
recently, the class-1,2 phases as well as other (conventional) 
phases of SU(3)-invariant spin chains were investigated from 
the SPT point of view.®® 

A remark is in order about the definition of the topolog¬ 
ical class. In contrast to the SU(2) case where all the ir¬ 
reducible representations are self-conjugate, we must distin¬ 
guish between an irreducible representation and its conjugate 
in S\J{N). The relation (18) suggests that if we have the edge 
state transforming under the projective representation TZ on 
the right edge, we necessarily have its conjugate TZ on the 
other. This means that when we talk about the topological 
class we must first fix which edge state we use to label the 
topological phases. Throughout this paper, we define the topo¬ 
logical class by the right edge state [i.e., Ug acting from the 
right in Eq. (18)]. Now it is easy to see that the SU(Af) VBS 
state introduced in Sec. IIC belongs to class-Af/2. 


C. Entanglement spectrum 

Remarkably, the above-mentioned difference in the pro¬ 
jective representation Ug can be seen in the entanglement 
spectrum.'"^ In order to define the entanglement spectrum, we 
first divide the system into two subsystems A and B. Then, 
the entanglement spectrum 0)} is defined through the 

Schmidt decomposition of the ground state \ip) of the entire 
system: 

IV-) = ® (19) 

CK —1 

where {!<))„)} and {|</>®)} are orthonormal basis sets for the 
subsystems satisfying = ^a /3 and the number x 

of finite ^q(< 00) defines the Schmidt number. 

According to Ref. 19, the entanglement spectrum of a given 
system exhibits a structure quite similar to that of the (energy) 
spectrum of the physical edge state of the same system and 
might be useful in characterizing topological states of matter. 
In one dimension, the edge states are not dispersive and we ex¬ 
pect a discrete set of degenerate levels to appear in the entan¬ 
glement spectrum reflecting the physical gapless edge modes. 
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In fact, in accordance with the degeneracy in the entangle¬ 
ment spectrum, the projective representation Ug assumes a 
block-diagonal structure,^® where each block corresponds to 
an irreducible representation of S\J{N) compatible with the 
topological class. For instance, in a ground state in the class-2 
topological phase of SU(4), each entanglement level should 
exhibit the degeneracy corresponding to an SU(4) irreducible 
representation with ny = 2 (mod 4). In Table I, the Young di¬ 
agrams as well as their dimensions are listed for some typical 
irreducible representations compatible with the class-2 topo¬ 
logical phase [i.e., ny = 2 (mod 4)]. 


I. VBS point 

To investigate the topological phase protected by PSU(4) 
symmetry, we begin with the simplest case. The ground state 
of the SU(4) VBS Hamiltonian (11) can be given exactly in the 
form of an and its entanglement spectrum is readily 

obtained by rendering the MPS into the canonical form (for 
the expressions of the matrices, see Appendix A). 

Reflecting the existence of the 6-dimensional (physical) 
edge states (ntop = ny = 2), the only entanglement level 
indeed is 6-fold degenerate indicating the class-2 phase:^° 
= log 6 (a = 1,..., 6; X = 6). This is in perfect agree¬ 
ment with the above argument. 


2. Heisenberg point 

In order to check if the ground state of the SU(4) Heisen¬ 
berg chain (1) is in the class-2 topological phase, we cal¬ 
culated the entanglement spectrum with the infinite time- 
evolving block decimation (iTEBD) algorithm.^'’^^ which en¬ 
ables us to directly access the entanglement spectrum. 

The simulations were done using the MPS with the bond 
dimensions up to 150 and the spectrum obtained is shown 
in Fig. 2. The degrees of degeneracy seen in Fig. 2 are 
{6,64,6, 50} from the bottom to the top. Clearly, this pat¬ 
tern perfectly fits into the dimensions in Table I; the edge state 
transform under the four (self-conjugate) irreducible represen¬ 
tations shown in Fig. 2. All these have ny = 2 (mod 4) and, 
from the discussion in Sec. Ill B, this ground state is classified 
as the topological class 2. 

Here a remark is in order. As the bosonic SU(4) Heisen¬ 
berg model (1) is obtained as the effective Hamiltonian in the 
Mott phase of the fermionic model (3), one may suspect that 
the same degeneracy structure could have been obtained for 
the original fermion model as well. However, this is not nec¬ 
essarily the case. In fact, in models where both bosonic and 
fermionic modes coexist, the entanglement spectrum contains 
the contribution from the fermionic sector as well as that from 
the bosonic one, and some of the levels may not obey the de¬ 
generacy rule that is obtained for the purely bosonic models. 
This is the reason why we simulated the effective bosonic 
model (8). 


TABLE I. Typical Young diagrams with the number of boxes ny = 2 
(mod 4) and their dimensions in SU(4). 


riY 

Young diagram 

dimension 

2 

d 

6 


□□ 

10 


10 

50 

64 



6 


70 



286 



FIG. 2. (Color online) Entanglement spectrum of an infinite 
SU(4) Heisenberg chain calculated by iTEBD. The degeneracy 
{6, 64,6, 50} may be understood in terms of the SU(4) irreducible 
representations shown in the figure, (inset) Zoom-up of the lowest 
six-fold-degenerate entanglement level. 


3. Continuity between Heisenberg and VBS points 

In the previous sections, we have seen, by inspecting the en¬ 
tanglement spectra, that the original SU(4) Heisenberg model 
(1) and the solvable SU(4) VBS model (11) share the same 
topological properties in common. Next, we consider adia¬ 
batic connection between the Heisenberg point and the solv¬ 
able VBS point to show that they belong to the same unique 
phase in the sense that they are connected to each other with¬ 
out quantum phase transitions.^^ To connect the two Hamilto- 
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nians, we use the following one-parameter family of Hamilto¬ 
nians 

n{a) 

= I + a 

i ^ 

( 20 ) 

where a is an interpolating parameter changing from 0 
[Heisenberg point; Eq. (1)] to 1 [VBS point; Eq.(ll)]. We 
calculated the entanglement spectrum of the ground state of 
•H(a) for a = 0.0, 0.1, 0.3, 0.5, 0.7, 0.9, and 1.0 with iTEBD 
and the results are shown in Fig. 3. It is evident that the struc¬ 
ture of the entanglement spectrum (including the six-fold de¬ 
generacy in the lowest level) is preserved all the way from the 
Heisenberg point up to the VBS point showing that the two 
models indeed belong to the same class-2 topological phase. 




Heisenberg 
14 


VBS 

4 


i 1 





1 

1 1 




1 

1 


. . .L . . 



1 

1 

■"Te "T" 




1 

1 

1 

1 

1 1 

1 1 




1 

1 

1 

1 

64 

1 1 

1 1 

1 1 




1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 1 

1 1 

1 6 1 




1 

1 

1 

1 

1 

1 

II 1 1 1 II 

II 1 1 1 II 


FIG. 3. (Color online) Evolution of entanglement spectrum as we 
interpolate between SU(4) Heisenberg model [(1); a = 0] and SU(4) 
VBS model [(11); a = 1]. Numbers shown next to the levels are 
degrees of degeneracy. 


IV. NON-LOCAL STRING ORDER PARAMETERS 


physical edge states to external perturbations, see Ref. 80). 
However, these order parameters are written directly in terms 
of the projective representation Ug and are not accessible in 
experiments in spite of their use in numerical simulations. 
Therefore, for the purpose of the detection of SPT phases in 
experiments, the characterization with order parameters, that 
are written in terms of measurable quantities, is still useful. 
In this section, we introduce a set of non-local string order pa¬ 
rameters for our SU(iV) spin system to characterize the topo¬ 
logical phases. 


A. Zjv xZiv and SPT phases 

In Ref. 81, a set of generalized string order parameters 
based on the symmetry was introduced for generic 

Zjv X ZAT-invariant systems and its connection to the Zat x Zjv 
SPT phases was discussed. As PS\J(N) and ZatxZat 
have the same cohomology group^*’^' iJ^(PSU(W), U(l)) = 
H'^{ZnxZ]s[= Ztv in common, we may expect that 
we can characterize our topological phase by using these 
string order parameters. In order to adapt the string order 
parameters, that was introduced in Ref. 81 in the context of 
Zat xZAT-invariant systems, to our SU(A^) case, we have to 
first identify the two commuting Zats in S\J{N). 

The construction of a pair of ZatS itself does not rely on a 
particular choice of the irreducible representation. In fact, we 
do not need the explicit expressions of the generators which 
depend on the choice of the basis and representation; the com¬ 
mutation relations among the generators suffice for our pur¬ 
pose. The most convenient way is to use the Cartan-Weyl 
basis {Ha, Ea} that satisfy^^ 

[Ha,Hb]=0, [Ha,Ea] = {a)aEa , 

3 

[Ea,,E_a,] = Y.^Oi)aHa , Tf = nSab, (21) 

a—1 

{a,h=l,...,N -1) 


In Sec. Ill C, we have seen that the structure of the entan¬ 
glement spectrum helps us to identify the topological class 
of a given ground state provided that we have enough infor¬ 
mation on the protecting symmetry of the system in advance. 
However, in general, the degeneracy structure alone does not 
uniquely identify the topological class. For instance, the class- 
2 phase of PSU(4)-symmetric systems has doubly-degenerate 
entanglement levels (see Appendix C), that are reminiscent of 
the Haldane phase protected by Z 2 XZ 2 , although these two 
phases are essentially different as we will see in Sec. V. Fur¬ 
thermore, despite some recent proposals,it is not very 
straightforward to directly measure entanglement in experi¬ 
ments. In fact, what is more fundamental in identifying SPT 
phases is the projective representation Ug. Therefore, “order 
parameters” that have more direct access to Ug is desirable. 

Several order parameters for SPT phases, including a 
gauge-invariant product of UgS, were proposed recently^^ (for 
discussion of the detection of SPTs using the response of the 


where a denotes the — N roots of SU(A^) normalized as 
|a| = v/2 which are generated by the simple roots (i = 
1,2,3). The normalization k depends on the representation 
and set to 1 for the A^-dimensional fundamental representation 
□ [e.g., At = 16 for the 20-dimensional representation ^ of 
SU(4) considered here]. In the actual calculations, one may 
use, e.g., the generators and the weights given in Sec. 13.1 of 
Ref. 82 with due modification of the normalization. 

Now let us look for the operators Gq and Gp that generate 
the two Zats. Regardless of N, the first generator Gq, which 
is diagonal and plays the role of in SU(2), is given simply 
by 

N-1 

Gq=Y. ’ ( 22 ) 

where are the A^ — 1 Cartan generators and p is the Weyl 
vector of SU(A^). The generator Gq has the following simple 
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commutation relations with the simple roots a: 


[Gq,E±^]=±E^ , (23) 

which guarantee integer-spaced eigenvalues of Gq (for 
the fundamental representation TV, they are essentially 
1, 2, • • • , TV). With this, the first Zn is generated as 

Q = cn exp f j , (24) 


where the phase cn has been introduced so that Q satisfy 
qN _ ^ rpj^g expression of the other generator Gp depends 
on TV and, in the following, we will explicitly work it out for 

TV = 4. 

The first Z 4 -generator Q is defined in terms of the two com¬ 
muting SU(4) generators (the Cartan generators) as 

Q = e*T exp = ^ ^^ 5 ) 

Gq = 2TTi + iJ2 • 

The generator Gq satisfies Eq. (23). On the other hand, the 
second Z 4 is generated by 


P = exp ( i^Gp ) , P'^ = 1 


Gp — „ 'y ] Efji + „ ( y ] Gjai Ti/Q,j_i_Q,2+Q3 j 


(26) 


_ 1 -E 

cy I / ^ ^ — Gtj ^ — Ol \— OL 2—0 




The summation twelve non-zero roots a 

of SU(4). Here we do not give the explicit expressions of the 
generators which depend on a particular choice of the basis, 
since giving the commutation relations ( 21 ) suffices to define 
Z 4 XZ 4 (see Supplementary Material for the expressions in 
a particular basis set that are more convenient for the actual 
calculations). It is important to note that the two operators 
Q and P constructed here generate Z 4 XZ 4 (i.e., [Q, P] — 0) 
only when the number of boxes in the Young diagram is an 
integer multiple of 4. In other words, what we have defined 
is the Z 4 XZ 4 subgroup of PSU(4). This is reminiscent of that 
the two TT-rotations along the x and z axes generate Z 2 XZ 2 
only for SO(3) ~ SU( 2 )/Z 2 . In Appendix B, we present the 
expressions of Gp and Gq for other TVs. 

Having explicitly constructed a ZatxZjv subgroup of 
PSU(TV), we now consider how the existence of this subgroup 
leads to (TV — 1) SPT phases. Consider the two Z^v generators 
P and Q satisfying 

(Q)^ = {P)^ = 1 ■ (27) 


If one requires that QP = PQ hold when both sides act on 
the MPS in question, one obtains from Eq. (18) 

A{m) = {UQUpUlul)A{m){UpUQUluG) . (29) 

When the MPS in question is pure and canonical, this implies 

UpUQ=e-^^‘^^UQUp . (30) 

On the other hand, combining (Up)^ = 1 and 

UpUqUpUq = obtained above, we obtain another 

relation; 


C/qC/^-1 = UqUI = U^-\UpUqUIuI)Uq 
= e-^’^»^U^-^UQ . 

Using (30), the right-hand side may be rewritten as 

= {e-^‘^^^fU^-‘^UQUp 

Therefore, we arrive at the conclusion that is the 

phase:^' 


(31) 


(32) 

Zn 


g*4>QP ^ ^ ^".op = 0,1,..., TV - 1) . (33) 

To see that <I>qp is in fact given by (27r/TV)ptop[= 
(27r/TV)nY], we just note that UqUp = e^^UpllQ for the 
TV-dimensional representation □ and that other representations 
are constructed by tensoring □ ny times. Eq. (33) implies that 
the exchange phase between Up and Uq carries the informa¬ 
tion on the topological class Ptop- 


B. Definition 

Next, we define another set of operators Xp and Xq satis¬ 
fying the following commutation relations with P and Q in¬ 
troduced in the previous section 

Q^XqQ = ujXq , P^XqP = Xq 

QtXpQ = Xp , P^XpP = uj-^Xp (w = e*^) 

for any irreducible representations of SU(TV). Using the com¬ 
mutation relations (21), one sees that the operators Xq and 
Xp for TV = 4 can be expressed by the SU(4) generators as 

Xq = —^{E-a.1 + E-a2 + P-az + P-on) 

v2 (35) 

Xp = H^- iH^ , 


As has been shown above, we can explicitly construct P and 
Q using the generators of SU(TV). By carefully choosing the 
gauge, we can always make the corresponding projective rep¬ 
resentations Up and Uq satisfy 

Ui = l, (Up)^ = (Uq)^ = 1, 

Up. = {Upr, UQ. = {UQr (n = l,...,TV-l) . 

(28) 


where q ;4 = —ai — a 2 — as and the normalization has been 
chosen such that Tr (XqXq) = Tr (XpXp). Erom these op¬ 
erators, we define the following string operators: 

Vp{m,nG) = 0^(1)” • • • Q\i - 1)” [Xp{i)y (36a) 

VQ{m, n- i) = (Xq(*))™ P(* + 1)- ■ • • P{Lr • 


(36b) 
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Then, the string-order parameters (SOP) are (infinite-distance 
limits of) the two-point functions of these string operators: 

Oi{m,n) = lim {Vp{m,n;i)Vp{m,n', j)) 

\i-j\y'oD 

(37a) 

02{m,n)= lim {VQ{m,n-,i)vX{m,n; j)) 

{0 < m,n < N) . 

(37b) 


where Uq is the projective representation of Q and the transfer 
matrices are defined as 

a,b—l 

d 

E [A*{a)]^j[A{b)]^^p{a\{XprQ-\b) . 

a,6—1 

(39) 

The ) denotes the largest left (right) eigenvector 

of the following transfer matrix: 

d 

[TqLm = E imUp iamb) . (40) 

a,6—1 



The subscripts 1 and 2 refer to the SOP corresponding to the 
two commuting Zjv’s (associated with Q and P, respectively). 

It is important to note that when the model is realized in the 
cold-atom system (3), the SOP Oi{m,n) are expressed only 
in terms of the local fermion numbers Ua^i = iCga,i + 
Cea In f^ct, the expressious (37a) involves only the di¬ 

agonal generators {Ha} [see, e.g., Eqs. (25) and (35)] which, 
when second-quantized, can be written only with the local 
fermion densities na,i- This property is desirable in view of 
future detection of the non-local order with the site-resolved- 
imaging techniques. 

As is seen in (37b), the second SOP 02{m, n) contain the 
off-diagonal generators [see Eqs. Eqs. (26) and (35)] and are 
more complicated; in order to express them in terms of the 
fermions, we first second-quantize the (off-diagonal) genera¬ 


tors, e.g., as 


Ea,i Cg^_j(^a)/37Cg7,i 


'-e/3, 


ii^a)l3-, 


where the 4x4 matrices Sa are four-dimensional fundamen¬ 
tal representations of the generators (see Supplementary 
Material for the expressions). Acting on the states in (2), the 
second-quantized generators Ea reproduce the ones appear¬ 
ing in (37b). 

The merit of using the SOP is that they carry the informa¬ 
tion on the projective representation Up and Uq that deter¬ 
mine the topological class^^’^® (see Sec. IV A). To show this, 
we first note that the SOP decouple into the product of the 
boundary contributions: 


Odm,n) E « 

—o ^ 0',oc;p,p 

OL^P 


Q 


Q 




y^pi 


Oi,L(w,n)e)i^R(m,n) , 


Using the properties of the canonical MPS,®® we can show 
that the right boundary term C9i R(m, n) in Eq. (38) satisfies 
the following identity^"’’™’*' (see Eig. 4): 

C>iR(m,n) =a;-'(™+"”“-)C>iR(m,n) {I = . 

(41) 

That is, if ^ 1 for some I, C>i_R(TO,n) = 

C9i(m, n) = 0 solely by symmetry. A similar identity is ob¬ 
tained for 02{m, n) as well. Then, these idendity imply that 
when both Oi(m, n) and 02 (w, n) are non-zero, the topolog¬ 
ical index ntop necessarily satisfies 

Eor A^ = 4, we can use the set of Oi 2 (m, n) with 

(m, n) = (1,3) (class-1), (2,1) (class-2), (1,1) (class-3) 

(43) 

to distinguish between the three topological phases (as well as 
one trivial one). In the SU(4) class-2 phase we discuss here, 
we expect 


01 . 2 ( 2 , 1 ) 340 , 

Oi.2(l, 3) = 01 , 2 ( 1 , 1 ) = 0. 


(44) 


In fact, for the solvable SU(4) VBS state discussed in 
Sec. IIC, we have 


jo (TO,n) = (l,3) 

Oi, 2 (to, n) = <1 1 (TO,n) = (2,1) (45) 

[0 (TO,n) = (l,l), 

which clearly indicate the class-2 topological phase. 

In general, we need a set of 2{N — 1) SOPs Oi, 2 (m,ri) 
to identify the PSU(Af) topological phases. Note that the non¬ 
vanishing SOP is the sufficient condition for the corresponding 
topological class. In other words, even if the system is in the 
topological phase, the corresponding SOP might be zero for 
some other special reasons. 


( 38 ) 
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FIG. 4. (Color online) Boundary term Oi^R(m,n) carries 
the information on the exchange phase a;"*"!’ between Up and 
Uq [see Eq. (30)]. Here a trivial identity (Xp)'^ = 

(pt)i |(p)*(Xp)"‘(pty| (py (;. arbitrary) has been used. 


C. Reflection 

In contrast to the SU(2) case where the operators Xp = 
and Xq — are hermitian (see Appendix B), Oi 2 (m,n) 
are not invariant under reflection symmetry I (with respect to 
a site or a bond) for SU(A^) with Af > 3. In fact, reflection T 
takes them to 


Oi, 2 {m,n) ^ Oi^ 2 {m,N - n)* , (46) 


where Oi 2 here are defined as 

Oi{m, n) 


\ I 

02 {m, n) 

(47) 



The new order parameters Oi 2 {m, n) look similar to the orig¬ 
inal SOP C7i 2 (rn,n) but are different in the relative posi¬ 
tion between the string and the end points [see Eqs. (37a) 
and (37b)]. Now one can repeat the preceding argument 
[see Eq. (38) and Eig. 4] on the boundary terms to obtain 
exactly the same selection rule (42). Therefore, one sees 
that when both Oi{m,n) and 02 {m,n) are non-vanishing 
in a given ground state lip), its parity partner I\ijj) has finite 
Oi{m,N — n) and 02{rn,N — n), and hence is in another 
topological phase characterized by Oi^ 2 {rn, N — n). Eor in¬ 
stance, the SU(4) class-1 topological phase characterized by 
01 , 2 ( 1 , 3) is the parity partner of the class-3 phase character¬ 


ized by 01 , 2 ( 1 ,4 — 3) = 01 , 2 ( 1 ,1) (Eig. 6; see Supplemen¬ 
tary Material for the explicit demonstration). 


D. Numerical results 


To demonstrate the use of the SOP in detecting the SU(A^) 
topological phases, we plot the value of the SOP Oi{m,n) 
for the model (20) obtained using iTEBD. Note that by the 
SU(4)-symmetry, we do not need to calculate 02(m, n). That 
Oi (2,1) is non-vanishing for 7((a) from a = 0 to a = 1 gives 
a strong evidence of the class-2 topological phase. 


Heisenberg 



1.0 


0.8 

T—1 

0.6 

tM 
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■ • 


VBS 


0.2 


0.4 0.6 

a 


0.8 


1.0 


FIG. 5. (Color online) Plot of SOP for 'H(a). C7i(2,1) is non-zero 
between the Heisenberg point (a = 0) and the VBS point (a = 1) 
giving additional evidence for the topological nature. 


E. Non-local transformation 


Before concluding this section, we give a remark on the 
connection between the SOP and the non-local unitary trans¬ 
formation (generalized Kennedy-Tasaki transformation) elim¬ 
inating the entanglement of the SPT phase that was first intro¬ 
duced in Refs. 62 and 63 for the SO(3)-based spin chains (see 
also Refs. 85 and 86 for recent discussions in the context of 
disentangler). A straightforward generalization of the above 
non-local unitary transformation to the PSU(A^) case may be 
given by** 


Ukt = exp 



Y,Gp{k)GQ{j) 

k<j 


(48) 


Then, it is easy to see that the string operators defined in 
Eqs. (36a) and (36b) transform (up to phase) as 

U^rj,Vp{m,n]i)UKT = Vp{m,m + n;i) 

UlTVQ{m,n;i)UKT = VQ{m,m + n;i) . 

This and Eq. (43) imply that repeated applications of (7 kt 
take the system from one topological phase to another (see 
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Fig. 6 ). In particular, the class-1 and 3 phases can be reduced representations 0] at site i onto the physical states \mi)\ 
to conventional phase with (spontaneously-broken) local or¬ 
ders, while the class -2 is not. [-P’(™j)]a /3 = i®\P)i ■ (52) 


IZ4XZ4-SSB) 

OaihO) 

Wiz 

|Z 


Oa{l, 0 ) UkT 

4XZ4-SSB) 

reflection I 


I class-3) 


[/T 

I class-2) < I class-2) 

0„{2.1) ^kt 0„(2,3) 

• t reflection I 


The metric matrix TZ creates the SU(lV)-singlet out of the two 
fractional objects \a)i and \fi)i+i on the adjacent site as (see 

Fig. 1): 

Isinglet) = 'R.ap\a)i\l3)^+i . (53) 

Then, we can show that Ui is given by the matrix TZ'. 

n^A{mi)n = {np{mi))n 

= P{mi)TZ = {TZP{m,)f = 


FIG. 6. (Color online) Generalized Kennedy-Tasaki transformation 
17 kt and three PSU(4) SPT phases. The class-2 state are mapped 
onto the state of the same topological class hy Ukt- Note that both 
01 , 2 ( 2 ,1) and 01 , 2 ( 2 , 3 ) characterize the same class-2 phase [see 
Eq.’(42)]. 


V. SYMMETRY REDUCTION 


In SPT phases, the list of possible topological phases is 
closely tied to the symmetry we impose on the system, and 
a phase which is topological under a certain symmetry may 
not be so when we consider a lower symmetry. Although the 
protecting symmetry PSU(iV) is automatically (i.e., without 
fine tuning) guaranteed almost perfectly in alkaline-earth cold 
fermions,it would be interesting, from the theoretical 
point of view, to consider the fate of the topological phases 
when PSU(A^) gets reduced. 


where (j)j is 0 when both P{mi) and TZ are symmetric/anti¬ 
symmetric, and TT otherwise [in our case, P{mi) are symmet¬ 
ric by construction]. Therefore, in order to know if Ux is an¬ 
tisymmetric or not, we have only to know how the SU(iV)- 
singlet is constructed out of \a)i and |/3)i+i. 

The SU(Af)-singlet is written as the following fully- 
antisymmetrized product of N = 2n states in the fundamental 
representation N; 


I singlet) 

{ikdk} 


~ ^Ukf'A.ik} t Gii2---G I'^n) ''' kin) 

partition I {A} 


. {Jk} 


(55) 


A. Systems only with reflection symmetry 


As the states inside the braces transform like 


We begin with the case where the PSU(A^) symmetry is bro¬ 
ken down to reflection symmetry with respect to the middle of 
a bond {link-parity I). As is emphasized in Refs. 14 and 66 , 
symmetry operations (whether local or non-local) which keep 
a given state (which we assume is represented as an MPS) 
invariant are expressed in the form of Eq. (18): 

A(m,) 4 A{m,Y = e^’^^UlA{m^)Ui , (50) 

where Ux satisfies t/J = ±C/x. Depending on the sign ap¬ 
pearing on the right-hand, there are two classes for systems 
with link-parity X (topological when —1 and trivial if 4 - 1 ).'"^ 
Now let us determine the sign for the SU(A^) VBS state 
shown in Fig. 1 . To this end, we hrst note that the MPS matri¬ 
ces A{mi) is written as 


n 



(56) 


the symmetry of TZ is encoded in that of the coeffi¬ 
cient C’{ij,};{jfc}(= il)- From the antisymmetry of 
imediatetely sees 

^{ikYAjk} = (-1)" C'bGTiG = (-^r^{3k}A^k} ■ (57) 

Therefore, under the symmtry-lowering perturbation, the 
class-2 SPT phase crosses over to the topological Haldane 
phase (a trivial phase) when n = N/2 = odd (even) (see 
Fig. 7). 


A{m,)=TZP{mi) , (51) 


B. ZjvXZat '—^ Z2XZ2 


where P{mi) is the projection operators from the two frac- As has been seen in Sec. IV A, we may regard the N — 
tional objects \a)i and |/3)i [in our SU(4) case they are two 6 1 PSU(A^) topological phases as protected by the subgroup 
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ZatxZat (see also Appendix C). In that case, the following 
commutation relation determines the topological classes:*^ 

UpUq = (ntop = 0,l,...,iV-l) . (58) 

As the above ZjsixZn contains the Z^xZi subgroup gener¬ 
ated by and when N is even, we may consider the 
symmetry reduction ZpsxZ-^ i—>■ Z 2 XZ 2 . From the relation 

UqiUpf/^ = {Upf'^Ug 

= {-ir’-^{Up)^/^Ug , 

one can easily see that the projective representations of the 
two Z 2 generators satisfy 

iUg)^^^{Up)^^^ = (-l)5^”'“p(C/p)^/2(t7Q)^/2 ^ (60) 

It is known'^’"^^ that in the presence of Z 2 xZ 2 -symmetry, 
the phase is topologically non-trivial when the projective 
representations of the two Z 2 S are anti-commuting, i.e., 
(—1) = — 1. This is possible only when 

= 2(2fc-f 1) (fc G Z) and Pjop = odd . (61) 

Since our SU(A^) (N: even) topological phase corresponds 
to ritop = N/2, it remains topological (i.e., Haldane phase) 
even after the symmetry gets reduced down to Z 2 x Z 2 when 
N — 2,6,10,.... When = 0 (mod 4), on the other 
hand, the topological phases considered here (Ptop = N/2) 
smoothly cross over to trivial ones. In Fig. 7, we summarize 
the crossover predicted here. 

VI. CONCLUSION AND OUTLOOK 

The possibility of realizing SU(Ai) symmetry using 
alkaline-earth cold atoms provides a new arena for the 
symmetry-protected topological phases. In this paper, we 
have studied the topological properties of the ground state 
of the SU(A^) Heisenberg chain (1) with the “spins” (2) at 
each site, especially for A^ = 4. This model is interesting 
as it is expected to describe the Mott-insulating region of the 
two-orbital SU(Ai) Hubbard model (3). From the analysis of 
the ground state of the solvable VBS Hamiltonian (11), we 
have suspected that the ground state of (1) belongs to one of 
the three topological phases predicted for the SU(4)-invariant 
systems. To substantiate this, we have calculated the entan¬ 
glement spectrum of an inhnite-size system with iTEBD and 
found that the degeneracy structure is perfectly consistent with 
that expected for the topological class (called class-2 in the 
text). In order to establish the adiabatic continuity between 
the Heisenberg model and the solvable VBS model, we have 
considered a simple one-parameter deformation 77 ( 0 ) of the 
Hamiltonian. The entanglement spectrum preserves its de¬ 
generacy structure all the way between the two models thus 
establishing the continuity. 

Then, we have investigated how the entanglement spectrum 
changes when the protecting symmetry gets lowered. Specif¬ 
ically, we have considered the situations where the original 



FIG. 7. (Color online) Fate of SPT phases protected by PSU(2n) 
when the symmetry is reduced down to link-parity T (Sec. V A) or 
Z 2 XZ 2 (Sec. VB). All these phases are labelled according to the 
irreducible representation(s) under which edge states transform. The 
label “n” of the SPT classes stands for the number ny (see the text) 
corresponding to the projective representations. When only the link- 
parity or Z 2 x Z 2 is imposed, only a part of them remains topological 
(warped arrows). 


SU(Ai) symmetry (which is perfect in the alkaline-earth cold 
atoms) is reduced to (i) link-parity and (ii) Z 2 XZ 2 . In both 
cases, the stability of our topological phase depends on the 
value of N\ when N = 4n -I- 2 (i.e., N = 2,6,10,...), we 
expect a crossover from our SU(A^) topological phase to the 
Haldane phase. 

Although the entanglement spectrum gives useful insights 
about the nature of topological phases, it may not fully char¬ 
acterize it. In fact, what is more fundamental is, at least from 
the group-cohomology point of view, the projective represen¬ 
tations which is the mathematical representation of the phys¬ 
ical edge states. The non-local string-order parameter (SOP) 
is appealing since it contains the information of the projec¬ 
tive representation in a manner that may be accessible in ex¬ 
periments. We have numerically calculated the SOP Oi{2, 1) 
with iTEBD and observed that it stays finite in our topological 
phase. This gives another support to our claim that the ground 
state of the SU(4) Heisenberg model (1) is in the class-2 topo¬ 
logical phase. 

At least two interesting questions remain to be answered. 
One is about the quantum phase transition(s) out of the topo¬ 
logical phase discussed here. In fact, an SU(A^) dimer¬ 
ized phase (called “spin-Peierls”) is observed numerically in 
Ref. 21 next to (i.e., on the smaller-U side of) the SPT phase. 
As the inclusion of higher-order terms in t/U may be mim¬ 
icked by adding terms higher order in to (1), we 

may include an extra term that favors dimerization to study 
the topological-dimerized quantum phase transition. 

Another interesting problem would be the nature of the 
strong-coupling (Mott) phase of the model (3) with odd- 
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N. In this case, the orbital degree of freedom is not fully 
quenched and we obtain an effective Hamiltonian different 
from (1), where the S\J{N) “spin” are highly entangled with 
the orbital degree of freedom.*^ As the nature of the effective 
Hamiltonian, which is reminiscent of the Kugel-Khomskii- 
type model^* for manganese, is not understood, it would be 
interesting to investigate it by the strategy used here. 
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Appendix A: MPS matrices for SU(4) VBS state 

In this appendix, we give the matrices necessary for the 
MPS representation of the SU(4) VBS state in Sec. IIC. 
The MPS for the SU(4) VBS state is given by the following 
product of six-dimensional matrices A{m.i) = Ar(mi) = 
r(mi)A (we follow the notations used in Ref. 71): 

|VBS) 

= ^ A(mi)A(m2) • • • A(mL)|mi) ® |m2) ® ■■ ■® |mi) , 
{m,} 

(Al) 

where the summation is taken over all the weights rrii = 
TO?) of the 20-dimensional representation of SU(4) 
and A is a diagonal matrix with non-negative diagonal ele¬ 
ments. Throughout this paper, we assume inhnite-size sys¬ 
tems where the MPS is given by inhnite-product of matrices 
A(mi). For several reasons, it is convenient to use the canoni¬ 
cal form of the above MPS, where the transfer matrix satisfies 
certain conditions. One possible choice of the canonical MPS 
is®®) 



/I 

0 

0 

0 

0 

0\ 


0 

1 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 


0 

0 

0 

1 

0 

0 


0 

0 

0 

0 

1 

0 


\0 

0 

0 

0 

0 

1/ 


(A2) 
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r(o,o,o)A = y6 
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Note that the diagonal elements of A are related to the entanglement spectrum {^a} by [A]qq, = e 
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/O 1 0 0 0 0 \ 

0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
00000-1 
Vo 0 0 0 0 0 / 

(A3e) 

/O 0 0 0 0 o\ 
0000-10 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
Vo 0 0 0 0 o/ 

(A3f) 

/O 0 0 0 0 1 \ 
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0 0 0 0 0 0 
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0 0 0 0 0 0 
Vo 0 0 0 0 o/ 
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Appendix B: Construction of Zjv xZat 

In Sec. IV A, we have explicitly constructed the Z 4 XZ 4 
subgroup of PSU(4) using the generators of the latter. Below, 
we give the expressions of the xZjv generators in terms of 
PSU(iV) for AT = 2 and 3. 

Regardless of N, the first generator Gq is given simply by 

N-l 

GQ=J2if^kHk, (Bl) 

where are the — 1 Cartan generators and p is the Weyl 
vector of SU(A^). The generator Gq has the following simple 
commutation relations with the simple roots a: 

[Gq ,e^] = e^, [Gq , a;_„] = , (B2) 

which guarantee integer-spaced eigenvalues of Gq (for 
the fundamental representation TV, they are essentially 


1, 2, • • • , TV). With this, the first Zjv is generated as 

Q = cn exp , (B3) 

where the phase cn has been introduced so that Q satisfy 
= 1 . 

The expression of the other generator Gp depends on TV. 


For Z 2 XZ 2 , we recover the well-known results®^’®^ 

Gq=pH = {H = 725^ p = 1/72) (B4a) 

Gp = -^iT„-^iT_„ = -5". (B4b) 

The operators Xp and Xq satisfying (34) are obtained as 

Xp = S\ Xq = S^. (B4c) 

For Z 3 XZ 3 , we have 

Gq = piHi + P 2 H 2 = 72TTi, (p = (72,0)) (B5a) 

3 

= (B5b) 
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where ai 2 are the simple roots of SU(3) and is defined 
by aa = —ai — a^- The operators Xp and Xq satisfying 
Eq. (34) are given by 

Xp = Hi- iH 2 (B 6 a) 

{E-a^ + E-a2 + E-a^} (c«3 = —Oil — 02 ) ■ 

(B 6 b) 

Appendix C: PSU(iV) and ZnxZn 

Since PSU(7V) and ZnxZm[C PSU(iV)] share 
the same cohomology group _ff^(PSU(A^), U(l)) = 

H'^{Z]\{XZ]y,U{l)) = Ztv, a phase which is topologi¬ 

cal under PSU(-/V) may remain so even if we weakly break 
PSU(iV) down to ZmxZm- As we have seen in Sec. Ill, 
when the system has the full PSU(A^)-symmetry, the en¬ 
tanglement spectrum exhibits the degeneracy pattern that is 
compatible with SU(A^)-symmetry. That is, the degeneracy 
of each entanglement level should find the corresponding 
entry in TABLE I. Now let us consider how the reduction 
of the symmetry down to a subgroup Zp[XZ^ changes the 
entanglement spectrum. 

As the unitary matrices Up,Q assume block-diagonal forms 
reflecting the structure of the entanglement levels, the relation 
(30) holds for each block corresponding to the degenerate en¬ 
tanglement levels A 

Up{\)Up{\) = e*^«-C/Q(A)C/p(A) = e*^”“pC/Q(A)C/p(A) . 

(Cl) 

This restricts the degree of degeneracy D\ of each entangle¬ 
ment level.Calculating the determinant of both sides of the 


above equation, one obtains 

det({7p(A)[/Q(A)) =detC/p(A)det{7Q(A) 

= (e*^"‘“p)^Mett7p(A)dett7Q(A) , 

(C2) 

which immediately implies = 1. When N and 

Ptop are mutually co-prime, D\ should be integer multiple of 
N. Otherwise, D\ of each level may be smaller. In particular, 
the entanglement spectrum of the class-1 PST phase exhibits 
the A^-fold degenerate structure for any N{'> 2), which is 
consistent with the results of the explicit calculation^’ for the 
SU(W) VBS chain based on another representation.^’’^^ 

Eor = 4 (Z 4 XZ 4 ), there are three topological phases 
(i) class -1 (Ptop = l)j (ii) class -2 (Ptop = 2 ), and (iii) 
class-3 (Ptop = 3). In the class-1 and 3 phases, D\ = 0 
(mod 4), while any even integers are allowed for D\ in the 
class-2 phase. Therefore, the degeneracy pattern observed 
in Sec. IIIC in general may be modified when we relax the 
full PSU(4) symmetry down to Z 4 XZ 4 , although the system 
still stays in the same phase. Eor instance, the lowest six- 
fold-degenerate level might be split into, e.g., three two-fold- 
degenerate levels. 

Eor instance, we may add to the original Hamiltonian (1) 
the following ZmxZ jy-invariant perturbation [see Eq. (34)] 

VAr=5tv^|(lpW)'^+(xl,(*))'^| , (C3) 

which is a generalization of the well-known single-ion 
anisotropy D '^^e usual spin chains [note Xp{i) = 

Xl,{i) = foriV = 2]. 
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